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Electron irradiation has been used to introduce point defects in a controlled way in the CuC>2 
planes of underdoped and optimally doped YBCO. This technique allows us to perform very accurate 
measurements of T c and of the residual resistivity in a wide range of defect contents %d down to 
T c = 0. The T c decrease does not follow the variation expected from pair breaking theories. The 
evolutions of T c and of the transition width with Xd emphasize the importance of phase fluctuations, 
at least for the highly damaged regime. These results open new questions about the evolution of the 
defect induced T c depression over the phase diagram of the cuprates. 



The establishment of the condensed state in high-T c 
superconductors (HTSC) is not fully understood. The 
occurence of thepseudogap in the underdoped part of the 
phase diagram [jj has raised the question of the coexis- 
tence of competing order parameters, or of the occurence 
of preformed pairs, with an eventual condensation of 
the pairs in a coherent superconducting state at T c . It 
has been proposed that, in these low carrier concentra- 
tion systems, T c could be determined by the phase stiff- 
ness of the order parameter The phase fluctuations 
could explain the occurence of a direct transition between 
superconducting and insulating states. An experimental 
approach which has been used extensively has been to 
study the influence of the disorder on both the normal 
and superconducting properties. In plane impurity sub- 
stitutions induce T c depression 3 , modification of the 
superfluid density n s 0.Q and local depression of the or- 
der parameter 0|. Whether the T c decrease results from 
pair breaking effects of the d-wave order parameter 
from the inhomogeneity of the order parameter (so called 
" swiss cheese model" Q) or from a reduction of the phase 
stiffness is still highly debated. 

In order to bring more accurate information on these 
issues we have undertaken careful studies of the influ- 
ence of controlled disorder in single crystals of HTSC 
cuprates. This can be achieved by high energy electron 
irradiations performed at low temperatures which intro- 
duce only point defects such as Cu and O vacancies in 
the Cu02 planes [TJ ll2> liil ■ The fact that a unique sin- 
gle crystal is progressively damaged allows us to study 
the influence of disorder with an accuracy impossible to 
attain with chemical substitutions. This gives us the op- 
portunity to study precisely the properties of samples 
with highly reduced T c , as we can indeed control the ir- 
radiation fluence to reach T c = 0. The results reported 
in this paper display remarquably simple phenomenolog- 
ical variations. We observe that T c quite unexpectedly 
decreases linearly with defect content down to T c = 0, in 
both underdoped and optimally doped YBCO, at vari- 
ance with Abrikosov Gorkov based theories. The anal- 
ysis of the variation of T c and of the transition width 
ST C reveals that the resistivity at T c is the relevant pa- 



rameter which determines the T c depression when T c ap- 
proaches zero. This agrees with the expected influence 
of quantum phase fluctuations. We shall discuss whether 
this approach can describe as well the variation with hole 
doping of the initial decrease of T c with increasing disor- 
der. 

The sing le crystals of YBCO are similar to those stud- 
ied in 14]. In-plane resistivities were measured by the 
Van der Pauw method. Special care has been taken to 
put the electrical contacts on the edges of the samples 
to insure an homogeneous flow of the electrical current 
through the samples. In the case of YBCO7, two different 
crystals of the same batch ( #1 and#2) have been studied 
[l5j . The irradiations were performed with 2.5MeV elec- 
trons in the low temperature facility of the Van der Graaff 
accelerator at the LSI (Ecole Polytechnique, Palaiseau). 
The samples were immersed in liquid H2 and the elec- 
tron flux was limited to 10 14 e/cm 2 /s to avoid heating of 
the samples during irradiation. The sample thicknesses 
(~ 20/im) were much smaller than the penetration depth 
of the electrons, which ensured a homogeneous damage 
throughout the samples 16]. 

For all samples we have found that Matthiessen's rule 
is well verified at high temperature, as the high T parts 
of the p{T) curves shift parallel to each other. This is 
exemplified in fig.l, which displays the T dependences of 
the in-plane resistivity p ab for YBCO7 #2. This confirms 
that even for very high defect content (xd ~9% in the 
planes 0]) the hole doping is not significantly modified 
as was already shown before for low Xd [HI Urn ■ 

For low T c samples upturns of the resistivity are dis- 
closed at low T and increase with increasing Xd ■ These 
" InT" contributions have been analysed in the case of 
YBCO6.6 US as being associated with a combination of 
single impurity scattering and localization effects. The 
latter become dominant for defect contents for which su- 
perconductivity is fully suppressed. 

To compare our results with pair breaking theories, it is 
useful to study the variation of T c with defect content Xd- 
The T c values reported hereafter were measured at the 
middle point of the resistive superconducting transition 
and the error bars were determined from the 10%-90% 
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FIG. 1: The resistivity is plotted versus temperature for the 
single crystal YBCO7 #2 irradiated at low temperature by 
2.5 MeV electrons. Data taken either after annealing at 150K 
(curve 6) or 300K (curve 7) are shown to coincide. 

values of the extrapolated normal state resistivity. As we 
shall see later on, the sharpness of the resistive transition 
is a good indication that the homogeneity of the damage 
in the sample is ensured. 

The' best estimate for :v,j is a priori the irradiation fiu- 
ence |19). For low Xd, data were taken in the irradiation 
set-up without annealing the samples above 150K. We 
found that T c and the normal state resistivity p ab (mea- 
sured at 150K) vary linearly with the irradiation fluence, 
so that Ap ab — p ab (lbQK) — p^ b ire (150K) represents Xd- 
In many cases the samples had to be taken to room T 
between irradiation runs. Although part of the defects 
were annealed in such processes, the p (T) curves were 
found to superimpose to those obtained before annealing 
(e.g. curves 6 and 7 in fig.l). Therefore Ap ab remains a 
good estimate for Xd- However, in YBCO6.6, the vari- 
ation with irradiation fluence of Ap ab , measured with- 
out annealing above 150K, departs from linearity for 
Ap ab > lOOpQ.cm, although T c still decreases linearly 
with the electron fluence. In this case Ap ab has then 
been replaced by Ap* ab , its linear extrapolation with flu- 
ence. We have therefore reported in figure 2 the data 
for T c versus both Ap ab and Ap* ab for YBC0 6 . 6 . For 
YBCO7 the deviation from linearity only occured for 
Ap ab > 200/if2. cm. The corresponding correction, which 
did not exceed 10% has not been performed in fig. 2. 

The linear variation of T c with defect content, down to 
T c = 0, is the most striking feature of the data displayed 
in fig. 2. This result contrasts with that expected from 
the standard Abrikosov-Gorkov formula which for d-wave 
superconductors gives T c as 0, |2(j : 



-ln(-^) =*(«+ )-*() (1) 

J-cO ^ ^ 



FIG. 2: Decrease of the critical temperature AT C as a func- 
tion of the defect content Xd- The latter is proportional to 
Ap ab measured at 150K for two crystals of YBCO7. In the 
case of YBCO6.6 the data (full circles) for Ap ab has to be 
corrected into Ap* ab (full diamonds) which is a best estimate 
of Xd as explained in the text. The AG curve (dotted line) is 
represented in the top part. The best linear fits (full lines) are 
also diplayed together with those for EK shown in fig. 4 using 
an analytical fit for p(T c ) versus Ap ab or Ap* ab (dashed lines). 



where ^(x) is the digamma function , a — h/ (2TrkBT c T) 
is the pair-breaking parameter, and 1/t oc Xd is the 
scattering rate in the normal state . The scale of the T c 
reduction is fixed by the ratio of H/t and ksT c so that 
the impurity pair breaking becomes more and more ef- 
ficient when T c decreases. This yields the well known 
negative curvature of the AG curve shown in fig 2, which 
is obviously not observed in the present data |2l| . Some 
published data for impurity substitutions 22] have been 
fitted with EqU Within the limited accuracy on the im- 
purity content and substitution site they could be fitted 
as well with a linear variation. 

Let us also consider the width of the superconduct- 
ing transition, which should reflect the inhomogeneities 
of the sample. The initial transition width ST C ~ 0.7K 
increases roughly linearly with Xd but never exceeds 5K, 
and then decreases when T c approaches zero, as can be 
seen in fig. 3. Inhomogeneities due to the irradiation pro- 
cess, such as electron energy losses through the sample, 
or a divergence of the electron beam at the level of the 
sample, etc... result in a macroscopic distribution of Xd 
with the full width 5xd = axd- Within the AG theory 
the transition width 5T C should scale with the derivative 
of the AG function. A large increase of 5T C would then 
be expected when T c reaches zero, contrary to the ex- 
perimental observation. From the observed quasi linear 
variation of AT C with defect content one should expect 
a linear increase of ST C down to T c ~ 0. The fact that 
5T C goes through a maximum rather indicates that super- 
conductivity becomes less sensitive to the distribution of 
defect content. Such a possibility will be further consid- 
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ered in the discussion hereafter. Let us point out at this stage that the observation in un- 

derdopedpure compounds of a linear relation between T c 
and n s |2q|. has led some authors to consider that this 
relation applies as well for impure samples ptiL I27I . The 
present observation of a linear decrease of T c definitely 
contradicts this simple guess. 

In a totally different approach, Emery and Kivelson 
(EK) argue that in low n s superconductors, T c might 
be determined by phase fluctuations of the order parame- 
ter. The temperature of the classical phase ordering T™ ax 
which is proportional to n s /m* can be much lower than 
the mean field critical temperature and is therefore an 
upper bound on the true T c . In that case the influence of 
disorder is to increase quantum phase fluctuations. They 
propose |10| that their magnitude is determined by the 
value of p(T c ), and that superconductivty disappears for 
a critical value pg, so that 



FIG. 3: Width of the superconducting 5T C defined as shown 
in the insert as a function of Ap ab for the two crystals of 
YBa2Cu307. The maximum of 5T C corresponds to a T c of 
~ 40K. The variation expected from AG, EK and linear 
dependence of T c with Sxa/xd ~ 7% are also plotted (see 
text) 



ln(Tr*/T c )=p(T c )/p Q He/T c ) (2) 
in which e is the energy scale of pairing interactions. 



Although AG theory is not applicable when localiza- 
tion corrections to the conductivity occur, the experi- 
mental deviations are already evident for a range of T c 
values for which these effects are negligible. A quantita- 
tive discrepancy with AG theory has been also evidenced, 
as the initial decrease of T c is much slower than predicted 
(HIS El - F° r YBCO7 for instance, theoretical estimates 
of AT C /Ap ab range from 0.7to 1.2 K/pfl. cm, a factor 
2 to 4 higher than the experimental value found here 
(~ 0.35 K/pTt. cm). It has been shown that this differ- 
ence can be taken into account in the framework of the 
AG model for a d-wave order parameter if one assumes 
an anisotropic impurity scattering |23j |. With appropri- 
ate parameters, one can explain the initial decrease of T c 
but with a similar dependence with Xd as the AG curve, 
which does not solve the contradiction with our results. 

A major reason for a breakdown of the AG theory 
might result from the very short coherence lengths in the 
high T c cuprates, which does not allow to assume a uni- 
form gap averaged over the disorder. This might allow 
to explain as well the variation of the superfiuid density 
n s with Xd as proposed by Franz et al |24| . Such a pos- 
sibility is reinforced by the recent STM data Q which 
reveal that the large depression of the DOS peaks, and 
therefore of n s , found nearby Zn impurities disappears 
on a lengthscale comparable to the coherence length £. 
As these regions of depressed superconductivity overlap 
for large Xd, a simplistic all or nothing model ( so called 
" swiss-cheese" ) might explain the negative curvature of 
n s (xd) which has been noticed by various authors 0,0. 



FIG. 4: The decrease of T c as a function of p(T c ) as defined 
in the insert of fig. 3 is compared to Eq.2. The best fits are 
obtained for pg equal to 600 and 1380 /j,Q.cm respectively for 
YBa 2 Cu 3 07 and YBaaCugOe.e- 

The fact that we have access in our experiment to T c 
values near zero allows us to test this dependence very 
precisely in fig. 4, where the data for T c are plotted ver- 
sus p(T c ). Reasonable fits of our data can be obtained 
with a large range of e values. As suggested by EK, we 
have therefore chosen a realistic value e ~ 1200-/^, the 
antifcrromagnetic exchange energy in YBCO7. As can 
be seen in fig. 4, the data can be well fitted by Eq.2 with 
T e max = 103^ and p Q = 600pn.cm for YBCO7 (76.6K 
and 1380/if2. cm for YBCOg.6)- Whatever the value cho- 
sen for e we always found a value of pq roughly two times 
larger in YBCOg.6 than in YBCO7 as it is mainly given 
by the value of p(T c ) corresponding to T c = 0. 
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We can as well determine the expected evolution of 
the transition width 5T C assuming 5xd — axd- We have 
done this for YBCO7, taking the derivative of EqElwith 
an analytical fit of the data for p(T c ) versus Ap ah . The 
variation of ST C calculated with the parameters deduced 
from the fit of Fig. 4 reproduces rather well the trend of 
the experimental data in Fig. 3 with a = 0.07. This 
analysis can conclusively be done solely from the data 
near T c = and allows to explain jointly the variation 
of T c and of 6T C , while all the alternative possibilities 
examined so far did not. So this emphasizes the im- 
portance of quantum phase fluctuations and of p(T c ) in 
determining the actual T c in this highly damaged region. 
Let us point out that the values for pg taken per CuC>2 
sheet [Isj (p 2 Q D = lOfcfi/D and 25M7/Q for YBCO7 and 
YBCO6.6) are much smaller than those observed in ion 
irradiated thin films (ref.12 in [13 )■ This is probably re- 
lated to the fact that irradiation damage in thin films 
usually induces much larger increase of resistivity than in 
single crystals. However the large variation of pg with 
hole doping found here is not anticipated by EK. This 
might stem from the fact that the analysis of EK is a 
priori more appropriate to describe underdoped cuprates 
with low superfluid density. One expects a behaviour 
completely different for overdoped materials for which T c 
should be determined by the mean-field transition tem- 
perature and disorder should be mainly pair breaking. 

This leads us as well to consider that our fits are not 
a proof that the analysis of EK applies for the initial 
decrease of T c particularly in optimally doped YBCO7 
which is at the borderline between underdoped and over- 
doped behaviour. It is interesting at this stage to re- 
mind the results obtained for the initial decrease of T c 
in cuprates with different hole dopings n^. We have 
shown that AT C /Ap 2 f increases steadily with nh [l4]|. 
Obviously the initial slope of Eq[3 that is AT C /Ap = 
— (T c //3g)ln(e/T c ) cannot explain this behaviour unless 



one introduces a large non physical reduction of pq with 
increasing hole content. If the number of carriers is taken 
as the number of holes rih , the initial slope is such that 
AT C scales as m*/r all over the phase diagram [l4|, which 
indicates that pair breaking in the d-wave superconduct- 
ing state might still contribute to the initial decrease of 
T c . This possibility is favored by the STM observation, in 
optimally doped Bi-2212, of quasiparticle states around 
Zn impurities, as expected from strong scattering the- 
ories [3. A reasonable explanation of the data would 
then be the occurence of a crossover from pair break- 
ing to phase fluctuations regime for high disorder. Such 
a crossover should progressively shift towards increasing 
disorder with increasing doping. 

In conclusion the possibility to create very homoge- 
neous disorder and well controlled defect contents Xd has 
allowed us to study carefully the depression of T c down to 
T c = in underdoped and optimally doped YBCO. We 
have shown that T c which decreases quasi linearly with 
Xd, does not scale with the reported variation of n s .The 
observation of the narrow transition width in YBCO7 for 
large Xd allowed us to suggest that p(T c ) is the relevant 
scattering rate which controls T c for high damage, as ex- 
pected from quantum phase fluctuations. The present 
experiments provide sufficiently detailed information to 
stimulate further theoretical calculations which should 
take into account both the amplitude and phase varia- 
tion of the order parameter. 

As for the highly damaged samples, we might anticipate 
that their superconducting properties should strongly dif- 
fer from those observed in weakly disordered compounds. 
Although this regime occurs near superconductor to insu- 
lator transition, further work is clearly needed to under- 
stand the relationship with the superconductor-insulator 
and the metal-insulator transition observed in other 2D 
structures. 
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